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Paths of tetragonal states between two phases of a material, such as bcc and fee, are called Bain paths. Two simple 
Bain paths can be defined in terms of special imposed stresses, one of which applies directly to strained epitaxial 
films. Each path goes far into the range of nonlinear elasticity and reaches a range of structural parameters in 
which the structure is inherently unstable. In this paper we identify and analyze the general properties of these 
paths by density functional theory. Special examples include vanadium, cobalt and copper, and the epitaxial path 
is used to identify an epitaxial film as related uniquely to a bulk phase. 



PACS numbers: 64.70.Kb, 61.50.Ks, 68.56.Eq, 62.20.Dc 

Pseudomorphic epitaxy of a cubic or tetragonal (001) 
film typically results in a strained tetragonal structure. 
If the stresses on the tetragonal state vanish and also 
the state corresponds to a local minimum of energy with 
respect to tetragonal deformations, the structure will be 
called a tetragonal phase. Such a phase will be stable or 
metastable depending on whether it has the lowest energy 
compared to other minima. Frequently metals have two 
tetragonal phases; sometimes both are cubic, e.g., bcc 
and foe Na and Rb sometimes one is cubic and the 
other phase is tetragonal, e.g., Ti and V ||||. They can, 
of course, also have phases with other structures, e.g., Ti 
also has a hep phase. 

Many paths can go from one tetragonal phase to the 
other. If the geometries along such a path have tetrag- 
onal symmetry and if they connect bcc and fee phases, 
the paths have been called Bain paths A purpose 
of the present work is to define and discuss a partic- 
ular Bain path which will be called the epitaxial Bain 
path (EBP). The EBP is produced by isotropic stress or 
strain in the (001) plane of tetragonal phases accompa- 
nied by vanishing stress perpendicular to the plane, such 
as pseudomorphic epitaxy produces on an (001) cubic 
or tetragonal film. Epitaxy provides a valuable means 
of stabilizing metastable phases and of putting phases 
under very large strains, both tensile and compressive in 
the plane of the epitaxy. The EBP of a material identifies 
the phase that has been strained, checks quantitatively 
the elastic behavior, which can be highly nonlinear, and 
predicts which phase of the material will form on a given 
substrate. Thus, in order to understand the properties of 
epitaxial films and new materials, the knowledge of the 
EBP is indispensable. 

A different Bain path has long been discussed, partic- 
ularly by Milstein Q , in which uniaxial stress is applied 
to a tetragonal state along the [001] axis accompanied by 
zero stress in the (001) plane; this path is conveniently 
called the uniaxial Bain path (UBP). We compare the 
two paths, EBP and UBP, which are both physically 
realizable. We show that both have the same lowest 



possible maximum energy or barrier energy of all Bain 
paths between the two tetragonal phases. However the 
EBP has a special value in relating strained tetragonal 
structures produced by epitaxy to particular tetragonal 
phases. Depending on where the tetragonal structure 
lies on the EBP, we shall show that a clear choice can be 
made which classifies the structure as a strained form of 
a particular phase. Thus for materials with both a bcc 
and a fee phase, the bulk (i.e., the interior) of an epitax- 
ial film can be uniquely classified as a strained bcc phase 
or a strained fee phase. Experimental geometry analyses 
of the bulk structure are frequently made by low-energy 
electron diffraction of ultrathin epitaxial films, which can 
then be placed on the EBP in the tetragonal plane to 
identify the phase of which it is the strained form. 

We employ density functional theory together with 
the full-potential linearized augmented plane wave 
method || . An important advantage of a first-principles 
calculation, i.e., a calculation with no adjustable parame- 
ters, is that it applies just as reliably to highly strained or 
even unstable geometries, as it does to the unconstrained 
ground state, i.e., to the stable phase. The same relia- 
bility cannot be asserted for empirically adjusted poten- 
tials, such as the empirical pseudopotential model used 
in Ref. 0, which is adjusted to reproduce ground-state 
properties of only certain structures. The first-principles 
calculation provides a check on the validity of empirical 
potentials far from the structural ground state and of 
linear elasticity theory. 

The calculation proceeds by evaluating the total energy 
E(a, c) in eV/atom as a function of the tetragonal lat- 
tice constants a and c. Contour lines of constant E are 
plotted in Fig.|l] for vanadium on the (a/abec - ^/K>cc) 
plane, where V = ca 2 /2 is the volume per atom. The 
contour lines show clearly two minima, one at the bcc 
structure and one at a tetragonal structure, and between 
them a saddle point at the fee structure. Thus vanadium 
has two fairly deep and well-separated minima and so is 
a good choice to illustrate the tetragonal paths between 
minima flq] . Since E has an extremum at the cubic struc- 
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tures HQ , the saddle point will occur exactly at the fee 
structure. 

Also plotted in Figjl] are the EBP and the UBP. The 
EBP is calculated from the minima of E(c) at constant 
a, which corresponds to the epitaxial situation, i.e., a is 
held fixed and c adjusts to minimize E(c) and thus makes 
the out-of-plane stress, i.e., the stress along [001], vanish. 
Conversely, the UBP calculation fixes c and lets a adjust 
to minimize E(a) making the in-plane stress vanish. Both 
paths go through the minima and the saddle point, hence 
the maximum energy on each path is the same. 

The saddle point is called in Ref. H "a special un- 
stressed tetragonal state ... at a local energy maxi- 
mum" [§). Obviously any path on the tetragonal plane 
that does not go through the saddle point must go 
through states of higher energy than the saddle point, 
since it will cross contours in either the upper or lower 
central sectors formed by the saddle-point contours at 
E = 0.29 eV/atom. Two common paths that do not go 
through the saddle point, hence have higher maximum 
energies than the EBP or UBP, are : 1. the constant 
volume paths at the volumes of the minima, which differ 
from the saddle-point volume. These paths have been 
called volume conserving Bain paths (|||; 2. the hy- 
drostatic pressure path, which has equal in-plane and 
out-of-plane stress components. 

Figure || plots E vs. a along the EBP and UBP. Also 
indicated is the range from a = 2.58 A (c/a = 1.6) to 
a = 2.66 A (c/a = 1.2) for which we find along the EPB 
the tetragonal structure to be inherently unstable J9| . In 
this range the slope of the EBP changes sign, and the 
stability condition 
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is violated. Here By are calculated from c\\ — 
(a 2 /V)d 2 E(a,c)/da 2 1 c 13 = {ac/V)d 2 E{a lC )/dadc, 
C33 = {c 2 /V)d 2 E(a,c)/dc for the strained lattice go|. 
Tetragonal states outside the range of Eq.(|]) cannot be 
stabilized with any applied stresses. When a path on the 
tetragonal plane reaches the boundaries of this unstable 
region, the tetragonal crystal becomes unstable and an 
epitaxial film will not be able to form many layers at the 
bulk spacing. Epitaxial films strained to a state near the 
unstable region should show large changes in the phonon 
spectra as elastic constants weaken, which would be in- 
teresting to study. 

One consequence of the existence of an unstable region 
is that a clear conceptual distinction can be made be- 
tween strained states of the phase at lower c/a (for vana- 
dium the bec phase) and strained states of the larger c/a 
phase (for vanadium at c/a = 1.83). If there was a con- 
tinuous path of allowed states between the two phases, 
there would not be a clear criterion for regarding an ob- 
served strained state as strained from a particular one of 
the two phases. 



A second consequence of the unstable region is that the 
strained system never gets to the saddle point, but trans- 
forms to another strained phase at an energy at or below 
the energy of the state where the path strikes the bound- 
ary of the unstable region. Thus for the EBP and UBP 
the breakdown energy is at least 0.05 eV/ atom lower than 
the saddle-point energy, and it is lower on the UBP. 

The EBP and UBP can be found for other metals by 
similar calculations; for Cu and Co our results are shown 
in Fig.|| and Fig.[|. The topology of the E(a, c) sur- 
face of cobalt is similar to the one of vanadium, with 
the two minima and the saddle point corresponding to 
the structures where the EBP and the UBP intersect. 
In this case, however, the stable cubic phase is face- 
centered (c/a = a/2) with of cc = 2.50 A, the bec structure 
corresponds to the saddle point (dbec = 2.83 A) and a 
metastable tetragonal structure is found at c/a = 0.92 
(a = 2.91 A) (TT). The energy barrier Et, cc ~ Ef cc is 
72 meV. 

The situation is different for Cu: the stable phase 
is at the fee structure (af cc = 2.53 A), but there is no 
other tetragonal metastable phase, and the energy sur- 
face E(a, c) shows a large plateau near the bec struc- 
ture for this reason, the EBP and the UBP share a 
common part for 1.13 < a/af cc < 1.15. 

Strained structures found by LEED are entered in the 
graphs; in each case the error limits intersect the EBP. 
Prediction of which phase will be produced and how 
strained it will be by a given substrate can be imme- 
diately made by entering the a of the substrate on the a 
axis of the EBP. 



For Co on Rh(001) [|T3| and on Cu(001) \L4\ the struc- 
tures lie on the left of the unstable region and thus both 
correspond to strained fee structures. However the error 
line for Co on Fe(001) [H is at a value of a above the 
unstable range and must be considered a strained form 
of the tetragonal phase of Co at c/a = 0.92; bec Co itself 
at c/a = 1 (a = 2.83 A) is unstable. 

The strained tetragonal structure for Cu on Pd(001) 
lies in the range of the EBP identified as unstable by 
our theory. We therefore disagree with the conclusions 
of a recent scanning tunneling microscopy (STM) ex- 
periment Eg , where such a film was "identified" as a 
metastable tetragonal phase. The structure measured 
in the high coverage regime (1 — 7 layers at T= 300 K, 
3 — 7 layers at higher temperatures) agreed within exper- 
imental error with the LEED analysis jl^j; however, the 
morphology of the islands observed in the STM pictures 
was taken as indication of a structure under small planar 
strain, and, therefore, of the existence of a Cu tetrago- 
nal phase with lattice constant close to apd- Our results 
predict that there is no tetragonal phase at c/a < t/2 
(and a > af cc ); we also show that the film is too far from 
the bec region of the EBP to be regarded as strained 
bec. Thus, we conclude that the properties of such 
Cu/Pd(001) film are not bulk like, but instead stabilized 
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by the (still ultrathin) film properties. 

In Fig.|| (as in Fig.|J) the epitaxial paths from linear 
elasticity theory are also drawn for fee Cu (and Co). 
Epitaxial strain on a cubic (001) surface in the linear 
elastic approximation says normal stress vanishes, hence 
8c/ c = —(2^/(1 — v))8a/a, where v is the Poisson ra- 
tio. A comparison with the calculated EBP in both fig- 
ures shows how the linear approximation fails to be valid 
for 8a/ a equals to few percent. In the analysis of the 
Cu/Pd(100) film, it is evident that the film is already far 
from the linear elastic regime. 

Generalizations of this calculation of physically real- 
izable paths between phases are possible with existing 
codes. Milstein JllJ has considered the instability of 
tetragonal structures with respect to orthorhombic defor- 
mations and found such instabilities on the UBP outside 
the range of c/a between the two minima; these insta- 
bilities probably exist on the EBP also, but have not as 
yet been studied. It looks practicable to explore with 
present codes even more general paths of homogeneous 
strain between tetragonal phases that make more gen- 
eral deformations of the unit mesh than orthorhombic 
deformations. Possibly a path of deformed lattices could 
be found with lower maximum energies, as the path of 
trigonal lattices of Ref. M . It would also be possible to 
find epitaxial paths of tetragonal states for slabs of a fi- 
nite number of layers, which would include relaxations of 
surface layers. 

In summary, paths of tetragonal states between two 
tetragonal phases include two particularly simple cases, 
namely, those produced by epitaxial (isotropic biaxial) 
stress and by uniaxial stress. Both cases correspond to 
maintaining certain surfaces stress-free and both have 
minimal maximum energies at a saddle point of energy 
that both pass through. Along these paths a range of in- 
stability exists, which separates states that can be clearly 
regarded as strained from one or the other of the two 
phases. The epitaxial path permits definite identifica- 
tion of the equilibrium phases of strained epitaxial films 
and could usefully be found for all metals with tetragonal 
phases. 
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FIG. 1. Contour lines of constant E in eV/atom for vana- 
dium on the tetragonal plane showing a bcc minimum (the 
zero of E) with at,cc = 2.93 A and Vbcc = 12.6 A (c/a 
a minimum at a = 2.41 A (c/a = 1.83) and V = 12.8 A 3 
a saddle point at a = 2.65 A (c/a = (2)5) and V = 13.1 
The epitaxial Bain path (EBP, full dots) and the uniaxial 
Bain path (UBP, open dots) are also shown. 
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FIG. 2. Plots of E in eV/atom along the Bain paths EBP 
and UBP for vanadium between the minima and through the 
saddle point. The unstable ranges of a are marked with arrows 
for both paths. 
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FIG. 3. Epitaxial Bain paths on the (a/a{ cc -V/V{ cc ) plane 
for Co, with LEED bulk structures for strained Co/Rh(001) 
@, Co/Cu(001) @, Co/Fe(001) (lBJ . The error limits 
are shown by a line through the point and come entirely 
from the estimated error in c from experimental analysis. 
The theoretical value of af cc for the fee phase is 2.50 A 
(Vfcc = 11.09 A ), which agrees with the experimental value 
p6| . The dot-dashed lines are obtained from linear elastic- 
ity theory, with elastic constants of the fee structure. The 
horizontal arrow mark the range along the EBP path where 
Eq.(p]) is violated. 
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FIG. 4. Same as previous figure, but for Cu. The LEED 
[0 (open triangle) and the STM Q (full triangle) bulk re- 
sults for Cu/Pd(001) are given as well. The theoretical value 
of afee for Cu is 2.53 A (Vfcc = 11.52 A ), while the exper- 
imental points use the experimental value of af cc = 2.56 A 
(Vfcc = 11.76 A 3 ). 
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